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In  this  paper  the  author  introduces  the  operator  :*  p  (y)A  with 

n 

higher  order  accuracy  for  approximation  to  the  differential  operator  D, 
where  ^  denotes  centered  difference  operator,  U  denotes  averaging 


operator , 


V>  -  I  Cm(u-I)".  c„-~  C  co  -  1  • 

m»0 

l 

A  class  of  new  many-knot  spline  basis  fi.  (P  (V))  N,  was  suggested.  The 

K,n  n  k 

smoothing  formulas 

f.  -  zr  /”  R.  (-7T-) f(t)dt  and  S.  f  -  I  f  fl 
k,n  h  ■'-«  k,n  h  k,n  *•  i  k,n 


are  discussed.  Accession  For 
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A  HEW  WAY  FOR  CONSTRUCTING  HIGHER  ORDER 
ACCURACY  SPLINE  SMOOTHING  FORMULAS 

Dong- Xu  pi* 

The  modern  mathematical  theory  of  spline  approximation  was  introduced  by 
I.  J.  Schoenberg  in  1946.  In  the  paper  [6]  he  studied  so-called  "B-spline 
basis".  A  B-spline  basis  can  be  normalized  in  various  ways.  One  of  them  is 
the  so  called  normalized  B-spline,  see  [2],  denoted  by  N^  k  for  the  B-spline 
function  of  degree  k  -  1  having  support  (x^,  xi+Jt) .  The  spline  smoothing 
formula  for  degree  k  -  1  to  an  arbitrary  function  f  can  be  represented  by 
S^f  -  £  f^N^  k*  This  approximation  has  been  used  in  curve  fitting 
successfully  [1].  [4]. 

In  order  to  improve  accuracy  of  the  smoothing  operator  S^,  the  author 

in  this  paper  suggests  a  new  spline  basis  denoted  ft  ,  instead  of  N.  .  . 

i,k,n 

Thus,  a  new  way  for  the  construction  of  spline  smoothing  formulas  is 

introduced.  I  prefer  calling  *  I  k  n  a  8mooth*n<T  operator  with 

grade  n  and  order  k.  In  here  when  n  “  k  Q  is  just  N^k  and  S,^0 

* 

is  the  same  as  Sk.  since  nf  e  +  <P^,  this  is  a  class  of  many-knot 
splines. 

Concerning  higher  order  accuracy  spline  smoothing  formulas, 

I.  J.  Schoenberg  [1946]  has  already  discussed  in  [6]  and  Z.  S.  Liang  studied 
the  many-knot  spline  smoothing  [4].  My  main  attempt  in  this  paper  is  to 
suggest  a  new  way  for  constructing  them. 
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1*  The  smoothing  operator 


Denote  the  centered  difference  operator  by  A  ,  defined  by 

n 

\f(x>  t-  f(x  +  |)  -  f(x  -  |)  . 

For  simplicity  let  h  ■  1,  and  A  »■  A  . 

The  B-spline  of  order  k  with  equally  spaced  knots  are  denoted  by 
and  it  can  be  represented  by 

\<x)  -  (AO*1)  6(x)  ,  (1.1) 

Ni,k<*)  »  0.2) 

where  D-1  is  the  integral  operator,  5  is  Dirac  6-function. 

It  is  our  purpose  to  find  a  more  exact  difference  approximation  to  the 
operator  D.  X  would  like  to  choose  following  ready-made  identity. 

Fact  1.1  ((5]  p.  43) 

.  .  2.  /  .  2  V  ...m2  (ml)  2m+1  . . 

1  +  y  )  .  /i  +  y  l  (-1)  ■— j— i  y  .  (1.3) 

ma0  ' 

Fact  1.2  The  following  expansion 


x 


sh  x  \ 


m«0 


C  <ch 
in 


D 


m 


holds.  Set 


then 


( 2m+1 ) ! I  ( 2m+ 1 ) ( 2m- 1)...3.1  , 

.  S  r  m  /  _ i  ***1  p  a  < 

2m+1  m-1  '  '  ( 2m+ 1 ) I J '  C0 


(1.4) 


2- 


<  ‘>!!'.Wnrrvm t»j»TSf.r': -s 


Proof  Prom  (1.3) 


.  7.  .  2k  ./,  .  2  r  .  ..m  ml  _m  2m+1 

log(y  +  /1+y)-/1  +  y  1  (-1)  77-77777  2  y 


( 2m+ 1)11 


'1  +  y2  I  c  2n 


x  x  /  2  /  2 

Let  y  *  sh  — ,  then  ch  —  *  »1  +y  ,  x  -  2  log(y  +  *1  +  y  ).  Thus 


x  ■  2  eh  "  I  Cm2”(shf) 
m«0 


-  2  ch  §  .h  f  I  C  (2  .h2  f) 
m«0  & 


sh  x  ][  C  (ch  x  -  1)”  . 

m«0  m 


Introduce  operators  E  and  defined  by 


E  f(x)  ;«  f (x  +  a)  , 


V(x)  2  (f(x  +  ?>  +  f(x  "  2))f  U  *’ 


U1  ' 


and  notice  the  relationships  between  those  operators  (see  [3],  p.  230) 

B  ■  e°,  ch  ^  .  y  , 

D  2  1  1 

.  .  D  2  ~  2  2  ~  2  -r 

2  sh  —  *  e  -  e  *  E  -  D  ■  A  . 

Use  ^  and  I  instead  of  x  and  1  in  (1.4) 
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D  ■  2  sh  r  I  C  ( ch  ~  -  I)' 

i  ra  i 

m*0 


-  I  C  (U-  I)"S 


n*0 


n 


i  C  (II  -  1)“2  +  R  , 
-  in  n 

m«0 


(1.5) 


where 


R  s-  2  sh  §  I  C  (ch  f  -  I)' 
n  z  w  2 

nr*n+1 


(1.6) 


Define  A^  as  the  first  part  of  (1.5),  i.e.. 


•r(n)  t  C  (P  -  I)”^  -  P  (P)A  , 

o  .  n  n 

m*0 


where 


P  <P) 
n 


l  C  (P  -  I) 
m=*=0 


m 


l  2~j  I  (-1>*"j<*>C  \  (J)E2 

j-0  m-j  3  i-0 


(1.7) 


In  the  general  case,  define 


zJn>  •  ■  VV\  • 


(1.8) 


This  is  a  collection  of  operators  approximate  to  D.  Beyond  doubt  Pn(1)  *  1# 
P0(P)  -  I. 

Fact  1.3  If  k  is  any  nonnegative  integer,  then  the  sum  of  all 

1  k 

coefficients  of  items  (Ph)  ln  the  expansion  equals  to  1. 

Notice  (1.6),  the  first  term  in  Rn  for  any  h 


(1.9) 


This  implies  the  following: 
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*-v- . ^’•■wawwa*--'?  r~'  ■ 


-I 


Theorem  1.1  Assume  that  f  e  C 


2n+3 


Then 


A^n>f(x)  -  Of <x)  -  2"3<n+1)Cn+1f(2n+3) 


.  t  .  .  n+1  ,  .  n+1  .  , 

where  5  e  [x - y  h,  x  +  -y  h] . 

Definition  We  call  the  operator  8 

grade  n  and  degree  f. 

It  is  to  be  noted  that  (A^^d  ^  is  just  as 

n 

Here  it  is  the  smoothing  operator  of  grade  0  and 


Fact  1.4  From  "nieorera  1.1,  if  g  @  P. 


2n+1 


on 


7<n)  -1 


n+1 


n+1 


A£“'d  g  -  g,  an  xe  [a+yh,  b-yh], 


(C)h2<n+1> 

smoothing  operator  with 

with  X.  J*  Schoenberg's, 
degree  k» 

[a,b],  then 


2.  A  class  of  many-knot  splines 


As  has  been  already  pointed  out,  the  B-spline  My  with  equally  spaced 
knots  (h  -  1)  is  the  result  of  the  Oth  grade  smoothing  operator  applied  to 
the  Dirac  6-function 

HJt  -  (AD“1)Mk_1  -  {^D-1>k6  .  (2.1) 

Mow  we  use  the  smoothing  operator  1  of  grade  n  for  the 

6-function  repeatedly.  He  can  define  a  class  of  spline  functions  which  as 
more  knots  than  Ny : 

\  _  «-  <P  (U))V  (2.2) 

K#n  n  K 

and 


‘W’  -  ■ 


k-1 


If  l  -  k,  then  n(x)  -  (A(n))k{^-t--->  which  has  knots 


5<n,K)  ,  .  cyi.ueit  j  .  0#1 


2(n+1 )k,  n  >  0  . 


We  often  take  (  ■  k  if  without  note. 

The  following  facts  can  be  proved  easily  in  the  same  way  as  the 
corresponding  facts  for  My. 

Fact  2.1: 


(D 

(2) 

(3) 

(4) 

(5) 


Vi.00  '  ak.n(-X'' 

fl.  (x)  »  0  for  all  |x|  >  , 

JCp  n  * 

D®a  (x)  -  (A(n))™n  (x),  0  <  m  <  k> 

k,n  k-ra,n 


-m 


Q  (x)  -  (A*n*  )*{xk+m”V(k+m-1)  l),  m  >  0; 


)  .  “k,n<X  +  «  - 

J—  * 


-6- 


-7- 


Then 


fk.„  '  i  C-  Vn<^’,(t,at  • 


(2.4) 


Theorem  2.3  If  f  e  C  (-“,“),  then 


r.  if^(x)!2dx  <  c  if(A)(x)i2dx 

"  Kjll 


(2.5) 


Proof  Take  the  derivative  of  order  t  for  (2.4 ),  and  the  integration  by 


parte,  and  notice  that 


|"--n-  -X'  P  (cos  x)|  <  1  • 

x  n 


If  f  is  a  discrete  valued  function  *  f(x^),  ■  xQ  +  ih,  then  a 


numerical  smoothing  formula  is  as  follows: 

-  l  •  ‘2'6> 

Formula  (2.6)  can  be  efficiently  applied  to  the  problems  of  curve  fitting  for 


discrete  data 


3 .  Examples 


In  this  section  some  discussions  which  are  helpful  for  applications  in 
practice  will  be  given. 

From  <2.2>,  with  *  -  k.  n  -  1,  k  -  1,2, 3,4.  we  show  the  particular 
representations  as  follows: 


|x|  <  ^  , 


ni,1<x>  - 


-  <  lx|  <  1  , 

lxl-1  , 


I X 1  >  1  ? 


n2,1<X) 


50  _  65  .  . 
36  36  |X|' 

«  _  49  - 

36  36  1 

22  15  ,  , 

3?  +  36  |X|' 


36  36 


|x|  <  i 


lx)  <  1  , 

1  <  lx)  <  |  , 

|  <  |x|  <  2  , 

I  x  |  >  2  » 


°3,1(X) 


462  878  2 

432  ”  432  X  ' 

858  1584  .  ,  706  2 

-  -rrr  +  777  x  ' 


432  432 
471  810 
432  ~  432 


,  ,  319  2 

lxl  +  432  X  ' 


627  +  —  Ixf 
432  432  1 

141  114  ... 


169  2 

432  X  ' 


lxl  <  ^  , 

\  <  lxl  <  1  , 

1  <  |x|  <  |  , 

I  <  lxl  <  2  , 


Ill  .  JLH  ixi  +  44-  x  ,  2  <  lxl  <  r  , 


ft  (x) 
4,1,x' 


-  - 


7920 

7776 

8444 

7776 

25212 

7776 

4152 

7776 

22440 

7776 

9060 

7776 

1308 

7776 

64 

7776 

0, 


20556 

2 

13059  .  . 

tttT  ,x| 

7776 

•  X 

+ 

3144 

|x| 

14268  2 

7776 

* 

7776  X 

53448 

|x| 

36036  2 

7776 

4 

7776  x 

11328 

|x| 

7956  2 

7776 

•f 

7776  X 

28560 

|x| 

11988  2 

7776 

7776  X 

9240 

|x| 

3132  2 

7776 

■4 

7776  X 

1128 

|x| 

324  2 

7776 

* 

7776  X 

48 

|v| 

12  2 

7776 

IXI 

4 

7776  X 

r  f  fxy 

7901  ,  ,3 
7776  *X*  ' 

•1--1  jxi3 
7776  ,X|  1 

J663  |x| 3 

7776  11  ' 

—  lx.3 

7776  ,X|  ' 


From  (2*2)  the  following  tables  are  given: 


Table  1: 


Ixl  <  ~  , 
\  <  lx|  <  1  , 

1  <  Ix|  <  |  , 

f  <  lx|  <  2  , 

2  <  |x|  <  |  , 

~  <  |x|  <  3  , 

3  *  I  x  |  <  ~  , 

\  <  lx|  <  4  , 

|x|  >  4  . 


10 


Assume  four  points  in  the  plane  are  given t 


S3, 

,1f 

S4,1f 

t 

x(t) 

y(t) 

I 

x(t) 

y(t) 

' 

1.6 

|  0.6S648  j 

1  i 

1.11870  i 

0.0804  | 

1.13213 

1.8 

\ 

!  0.89167 

1.08685 

0.83711 

1.09813 

2.0 

!  1.02315 

i 

1.02315 

1.00514 

1.00514 

2.2 

\ 

;  1.08685 

j  0.89167 

!  1.09813 

|  0.83711 

2.4 

!  1.11870 

|  0.65648 

1.13213 

|  0.61804 
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